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METHOD FOR M0DBU.IN6 THREE DIMENSIONAL OBJECTS 
AND SIMULATION OF FLUID FLOW 

This invention relates to a method for modelling solid 
5 objects, particularly for use in the sixmilation of fluid 
flow, to be used for example to simulate prototypes before 
production. In a preferred embodiment the method is used 
in the design of articles to be manufactured by injection 
molding, preferably from molten plastic materials. 

10 

The modelling of solid objects is employed in various 
fields. Such modelling is used, for example, in the 
simulation of injection molding. In that field, it is 
widely recognized that the filling and packing phases of 

15 injection molding have a significant effect on the visual 
and mechanical properties of a molded object. Simulation 
is employed to analyse proposed shapes and Injection 
points, and thus the final quality of the ultimate article. 
A requirement of any injection mold is that it can be 

20 filled with molten polymer given the pressure limits of a 
real injection molding machine. Simulation can provide 
information as to whether the mold can be filled and the 
fill pattern that will be achieved. By using simulation, it 
is possible to determine optimum gate locations and 

25 processing conditions. It is possible to predict the 

location of weld lines and air traps. Economic benefit is 
derived from simulation because problems can be predicted 
and solutions tested prior to the actual creation of the 
mold. This eliminates costly re-working and decreases the 

30 time required to get an object into production. 

Simulation technology has been developed ainri generally uses 
finite element /finite difference techniques to solve the 
governing equations of fluid flow and heat transfer, in 
35 order to minimize the time required for analysis and hence 
the required computer resources, the Bele-Shaw 
approximation is invoked to simplify the governing 



wo 98/43179 



- 2 - 



PCr/AU9S/00130 



equations. It has been £ound that this siiiq;>lif ication 
provides sufficient accuracy for injection xoolding but does 
create the need for specific modelling of the coi^putational 
domain. 

5 

Injection molding is an excellent process for repetitively 
manufacturing large numbers of objects or parts having 
complicated geometries. A characteristic of injection 
molded components is that the thickness of the wall is 
10 generally a small fraction of the overall length of the 
component. In view of the low thermal conductivity of 
plastics, this physical characteristic is essential to 
achieve the rapid cycle times that make the process so 
attractive. 

15 

The flow of melt in an injection mold is determined by the 
familiar conservation laws of fluid mechanics. Solution of 
the equations in their full generality presents several 
practical problems. Owing to the characteristically thin 

20 walls of molded components, however, it is possible to make 
some reasonable assumptions that lead to a simplification 
of the governing equations. These simplified equations 
describe v^t is called Hele-Shaw flow and may be readily 
solved in complex geometries using an appropriate numerical 

25 technique such as the finite element and/ or finite 
difference method. 

Injection molding simulation is now routinely regarded as a 
desirable aspect of plastic part design. Similarly, 

30 improved cosiputer aided drafting (CAD) technology has led 
to the widespread use of surface and solid imodelling. 
Associated advantages of this are the ability to better 
visualize an object, to. use numerical cutting, and the 
possibility of achieving more concurrency in engineering 

35 design and manufacture. When using the Hele-Shaw 

approximation, plastic CAE analysis still requires the use 
of a surface model, representing the midplane of the real 
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component, which is then meshed with triangular or 
quadrilateral elements to which suitable thicknesses are 
ascribed. The ^preparation of such a mesh can take a 
considerable amount o£ time, and requires sxibstantial user 
5 input; owing to the labour intensive nature of this step, 
model preparation requires the greatest share of time when 
performing a molding simulation and makes this technique 
time consuming, in addition, as model preparation is an 
interactive task, it has a higher cost associated with it 
10 than simply running a. computer program. 

This conventional technique for the simulation of three 

dimensional bodies is illustrated in figures 1 to 3 of the 

accompanying drawings in which: 
15 Figure 1 is a representative exaii^le of an 

article to be manufactured by injection molding in the form 

of a T-shaped component; 

Figure 2 is a midplane representation of the 

component of figure 1 that is meshed; and 
20 Figure 3 is a view of the major results from a 

conventional analysis showing the position of the advancing 

melt at given times. 

One solution to the above shortcomings is to avoid the use 
25 of the Hele-Shaw equations luid solve the governing 

equations in their full generality. This has inherent 
problems owing to the thin walled nature of injection 
molded objects and parts. To perform such an analysis, the 
region representing the miold cavity into v^ich molten 
30 polymer will be injected must be divided into small sub- 
domains called elements • Usually these elements €ure of 
tetrahedral or hexahedral shape. This process of sub- 
division is called meshing and the resultant network of 
tetrahedra or hexahedra the mesh. Owing to the cosq^licated 
35 shape of many injection molded objects and parts it is 
generally not possible to automatically mesh the cavity 
with hexahedral elements. It is possible, however, to mesh 



wo 98/43179 



- 4 - 



PCT/AU98/0ai30 



the domain automatically with tetrahedral elements. The 
thin walled nature of injection molded objects parts 
means that the plastic is subject to a huge thermal 
gradient in the thickness direction of the component. This 
5 requires that there be a reasonable number of elements 

through the thickness. Using existing meshing technology, 
the result is a mesh consisting of hundreds of thousands or 
even millions of elements. The high number of elements 
makes the problem Intractable for any but the fastest super 
10 computers. These are rarely found in Industry/ being 

extremely costly to purchase and maintain. Thus, although 
three dimensional simculation provides a solution that 
avoids the requirement of a midplane model, it is not as 
yet a practical solution. 

15 

It is an object of the present invention, therefore, to 
provide a method for the simulation of flow in a three 
dimensional object that can produce simulations 
substantially automatically, without requiring the solution 
20 of the governing equations in their full generality. 

According to a first broad aspect of the present invention, 
therefore, there is provided a method for simulating fluid 
flow within a three dimensional object having first ^^n^ 

25 second generally opposed surfaces Including; 

matching each element of said first surface with 
an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 

3d matched elexnents . and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 

35 all locations on said first and second surfaces from which 
flow may eman ate such that resulting flow fronts along said 
first and second surfaces are synchronized. 
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It should be noted that, in this context, generally opposed 
surfaces may be parallel, inclined to each other, meet at 
an acute angle or otherwise, and need not be planar, 
5 . provided that a reasonable thickness or thicknesses may be 
assigned to the space between the surfaces. 

Preferably said injection point is first linked to all said 
locations from substantially the commencement of said flow 
10 analysis. 

Preferably said injection point remains so linked at 
substantially all times in said flow analysis subsequent to 
said first linking. 

15 

Thus, although. the injection point must be linked to all 
these locations at some point in the flow analysis, a^*^ 
essentially continuously thereafter, the injection point 
may remain unlinked initially* Although this may decrease 
20 the accuracy of the final analysis adversely, it may 

nevertheless allow the analysis to yield useful results. 

Preferably said injection is one of a plurality of 
injection points. 

25 

Preferably said synchronization of said flow fronts is 
checked periodically. 

Preferably said flow fronts are synchronized if found not 
30 to be or no longer synchronized. 

Preferably said checking is performed at each time step. 

Thus, although for simple objects it is sufficient to 
35 s3^chronize flow fronts by linking the injection point or 
points to all locations, more complex parts require this 
synchronization at each calculational time step. 
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Preferably said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 
surfaces . 

5 

Thus, where the existing techni<xae8 utili2Ee a single, 
midplane representation of the object in which flow is to 
be modelled, and perform a simulation with this 
representation, the method according to the present 

10 invention utilizes only the outer surfaces defining the 

three dimensional object to create a computational domain* 
These correspond to representations (in a preferred 
embodiment, meshed representations) of the domain in which 
flow is to be simculated, azid would comprise, for exaaqple, 

15 meshed representations of the top and bottom surfaces of a 
part* Thus, in such an example the invention could be said 
to utilize an outer s)cin mesh rather than a midplane mesh. 
Elements of the two surfaces are matched, based on the 
ability to identify a thickness between such elements * An 

20 analysis, substantially along conventional lines (by meazis, 
for exai^ple, of the Hele-Shaw equations), is then performed 
of the flow in each of these domains in which flow is to be 
simulated, but linked to ensure fidelity with the physical 
reality being modelled. 

25 

PreferaO^ly axxy unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 
of adjacent matched elements where such adjacent matched 
30 elements exist, or of adjacent unmatched elements where 
such adjacent matched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

Thus, eventually all matched and unmatched elements of the 
35 first and second surfaces may be assigned thicknesses. 

Preferably each element of an edge surface, being a surface 
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between said first and second surfaces, and adjacent to 
either of said first or second surface is assigned a 
thickness proportional to the thickness of the element of 
said first or second surface to which said element of said 
5 edge surface is adjacent. 

Preferably each said element of an edge surface is assigned 
a thickness between 0.5 and 1.5 times said thickness of the 
element of said first or second surface to which said 
10 element of said edge surface is adjacent. 

Preferably each said element of an edge surface is assigned 
a thickness between 0.7 and 0.9 times said thickness of the 
element of said first or second surface to which said 
15 element of said edge surface is adjacent. 

Still more preferably each said element of an edge surface 
is assigned a thickness 0.75 times said thickness of the 
element of said first or second surface to which said 
20 element of said edge surface is adjacent. 

Preferably each element of an edge surface not adjacent to 
said first or second surface is assigned a thiclmess being 
the average of the thicknesses of adjacent elements of said 
25 edge surface that have been assigned thicknesses. 

Assigning a thickness, therefore, to every element of both 
surfaces and of edge surfaces ensures that the greatest 
accuracy will be achieved in the simulation. 

30 

Preferably flow is simulated at a rate directly 
proportional to a desired flow rate for the object. 

The desired flow rate will usually be a flow rate selected 
35 by a user based simply on the volume of the object to be 
filled and the desired fill-time. 
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Preferably said rate is proportiostal to tlxe ratio of 
computational domain volume of said object to real volume 
of said object. 

5 More preferably said rate is substantially equal to the 
ratio of said computational domain volume to said real 
volume. 

This approach may be used to compensate for the 
lb modification of the effective volume of the object being 
modelled, owing to the use of two domains (for example, 
associated with the top and bottom surfaces ) , rather than 
one meshed domain as in the midplane approach. The 
computational domain volume is thus this effective volume 
15 used for the purposes of modelling the object, which most 
preferably will be approximately double the real volume. 

Preferably the method is performed with first and second 
representations of said first and second surfaces 
20 respectively comprising first and second meshes or lattices 
respectively, wherein said elements are interstices of said 
first and second meshes or lattices. 

Preferably said elements are triangular or quadrilateral. 

25 

Preferably, when said elements are triangular, said 
. elements are substantially equilateral. 

Preferably said method includes creating said first 
30 secoxid representations. 

In one embodiment the method includes creating improved 
representations of said first and second surfaces, \^ereby 
said el emen ts are elements of said improved r^resentations 
35 and said method is performed with said improved 
representations . 
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Preferably said first and second representations are, or 
are a part of, a representation or representations for 
stereolithography of said object. 

5 Thus, computer representations of the two surfaces may be 
provided as input to the method, may be created by the 
method, or - if provided with representations of the 
surfaces - the method may create improved representations . 
As described above, particularly preferred representations 
10 include those with elements comprising equilateral 
triangles . 

In one embodiment of the invention the method described 
above is performed by a computer ruxming a computer program 
15 encoding said method for simulating fluid flow. 



Preferably said method optionally includes corrections for 
non- isothermal temperature fields and/or non-Newtonian 
fluids. 

20 

Thus, known techniques for including the effects of non- 
isothermal temperature fields and non-Newtonian fluids may 
be included so that more physically faithful simulations 
may be performed, as well as faster simulations without 
25 these corrections when speed is desirable, even at the 
expense of precision. 

According to a second broad aspect of the present invention 

there is provided a method for simulating fluid flow within 
30 a three dimensional object having first and second 

generally opposed surfaces including: 

providing or creating first and second 

representations of said first and second surfaces 

respectively, 
35 creating first and second improved 

representations from said first and second representations 

respectively. 



wo 98/43179 



- 10 



PCT/AV98/00130 



matching each element of said first improved 
representation of said first surface with an element of 
said second improved representation of said second surface 
between which a reasonable thickness may be dief ined/ 
5 wherein matched elements of said first improved 

representation constitute a first set of matched elements 
and matched elements of said second improved representation 
constitute a second set of matched elements, 
specifying a fluid injection point, 

10 performing a flow analysis using each set of said 

matched elements, whereby said injection point is linked to 
all locations on said first and second improved 
representations from which flow may emanate such that 
resulting flow fronts along said first and second improved 

15 representations are synchronized* 

Preferably the first and second representations are, or are 
a part of, a representation or representatiozis for 
stereolithography of said object. 

20 

Preferably the first and secoxui improved representations 
comprise small equilateral triangular elements. 

According to a third broad aspect of the invention there is 
25 provided a computer provided with or running a coiQputer 
program encoding the method for simulatix^r fluid flow 
described above. 

In a fourth broad aspect of the invention there is provided 
30 a computer storage medium provided with a computer program 
embodying the method for simulating fluid flow described 
above. 

According to a fifth broad aspect of the present invention 
35 there is provided a method for modelling a three 
dimensional object including: 

specifying first and second generally opposed 
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surfaces of said object, 

forming first and second representations of said 
first and second surfaces respectively, wherein said first 
and second representations each conprise a plurality of 
5 elements, 

matching pairs of elements of said first and 
second surfaces between which a reasonable thickness may be 
defined. 

10 Preferably the first and second representations coi^prise 
first and second meshes or lattices respectively, wherein 
said elements are interstices of said first and second 
meshes or lattices. 

15 Thus, the elements of the first representation are 

interstices of the first mesh, and the elements of the 
second representation are interstices of the second mesh. 

Preferably said elements are triangular or quadrilateral. 

20 

Preferably, when said elements are tricuigular, said 
elements are substantially eczuilateral. 

Preferably each element of each of said matched pairs of 
25 elements is assigned respectively said thickness. 

Preferably any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 
30 of adjacent matched elements where such adjacent matched 
elements exist, or of adjacent unmatched elements where 
such adjacent matched elements do not exist and said 
adjacent unmatched elements have been assigned thicknesses. 

35 Preferably each element of an edge surface, being a surface 
between said first cmd second surfaces, and adjacent to 
either of said first or second surface is assigned a 
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thickness proportional to the thickness of the eleaent of 
said first or second surface to which said element of said 
edge surface is adjacent. 

5 Preferably each said element of an edge surface is assigned 
a thickness between 0.5 and 1.5 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

10 Preferably each said element of an edge surface is assigned 
a thickness between 0.7 and 0.9 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

15 Still more preferably each said element of an edge surface 
is assigned a thickness 0.75 times said thickness of the 
el emen t of said first or second surface to which said 
el emen t of said edge surface is adjacent. 

20 Preferably each element of an edge surface not adjacent to 
said first or second surface is assigned a thickness being 
the average of the thicknesses of adjacent elements of said 
edge surface that have been assigned thicknesses. 

25 According to a sixth broad aspect of the present invention 
there is provided a method for sixmilating fluid flow within 
a three dimensional object having first and second 
generally opposed surfaces including: 

matching each element of said first surface with 

30 an element of said second surface between which a 

reasonable thickness nay be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements, 

35 specifying a fluid injection point, 

performing a flow analysis using each set of said 
matched elements, and 
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synchronizing flow fronts resulting from said 
flow analysis along said first and second surfaces. 

Preferably said flow fronts are synchronized from 
5 substantially the commencement of said flow analysis. 

Alternatively said flow fronts are first synchronized after 
the commencement of said flow analysis. 

10 According to a specific aspect of the present invention 

there is provided a method for modelling three dimcmsional 
objects by finite element azialysis including: 

modifying a representation for stereolithography 
of said object to a representation suitable for finite 

15 element analysis. 

In one embodiment the method includes creating inqproved 
representations of said first and second surfaces, whereby 
said elem en ts are elements of said improved representations 
20 and said method is performed with said improved 
representations . 

According to a further broad aspect of the present 
invention there is provided a method for simulating fluid 
25 flow within a three dimensional object having first ^r\^ 
second generally opposed surfaces including; 

matching each element of said first surface with 
an el emen t of said second surf ace between which a 
reasonable thickness may be defined, wherein matched 
30 elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using said first set 
35 of matched elements, 

adapting said flow analysis to said second set of 
matched elements, and 
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synchronizing flow fronts resulting from said 
flow analysis and said adaptation of said flow analysis 
along said first and second surfaces* 

5 TfauB, it is possible in some cases to perform the method 
using only one of the sets of elements, and mapping the 
result onto the other set, and performing minor adaptation 
of the analysis to allow for any differences between the 
two surfaces. 

10 

Preferably the method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 
respectively, wherein said elements are interstices of said 
15 first and second meshes or lattices. 

PrefercO^ly the elem e nts are tricuigular, and most preferably 
equilateral triangles. 

20 Prefereibly the elements are quadrilateral. 

Preferably the method includes creating said first and 
second representations, and more preferably creating 
improved represehtatiozis, either from the original surfaces 
25 or from the first and second representations, and 

performing the method with elements of the representations 
or Improved representations. 

In all the above methods according to the present invention 
30 aspects it is preferred that the synchronization con5>rises 
matching pressure and temperature. 

According to particular aspect of the present invention 
there is provided a method for simulating fluid flow within 
35 a three dimexisional object having first and second 
generally opposed surfaces including: 

matching each element of said first surface with 
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an el emen t of said secozid surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
5 surface constitute a second set of matched elements, 

specifying a plurality o£ fluid injection points, 
performing a flow analysis using each set of said 
matched elements, vAereby said injection points are linked 
to all locations on said first and second surfaces from 
10 which flow may emanate such that resulting flow fronts 
along said first and second surfaces have matched flow 
rates. 

A preferred embodiment of the invention will be described, 
15 by way of example, with reference to the accoiqpanying 
drawings, in which: 

Figure 4 is an example of a stereolithography 
mesh of the T-shaped object of figure 1; 

Figure 5a shows an initial mesh used in a 
20 remeshing method according to the present invention; 

Figures 5b-f show the mesh of figure 5a 
progressively remeshed according to the remeshing method; 

Figure 6 shows the mesh of figure 4, remeshed 
according to the preferred embodiment of the present 
25 invention; 

Figure 7a is a view of a mesh produced by 
conventional modelling techniques for a flat plate; 

Figure 7b is a view of a corresponding mesh 
produced for a flat plate according to a preferred 
30 embodiment of the present invention; 

Figure 8 illustrates surface matching for a flat 
plate according to the preferred embodiment; 

Figure 9 is a view similar to figure 8 
illustrating in cross-section the edges, matched surfaces 
35 and tinmatched surfaces of another T-shaped object according 
to the preferred embodiment; 

Figure 10 is a cross -sectional view Illustrating 
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the edges, matched surfaces, matched tapered surfaces, 
unmatched surfaces and unmatched edges of yet another 
shaped object according to the preferred embodiment; 

Figure 11 is a view of a simple object showing 
5 injection points; 

Figure 12 is a view of an object comprising 
intersecting plates with a user selected injection point 
according to the preferred embodiment; 

Figure 13 is a enlarged view of the object of 
10 figure 12 showing the injection points required to allow 

flow to emanate correctly from the selected point according 
to the preferred einbodiment ; 

Figure 14 is a cross-section of a plate, with a 
flow front advancing from the left and leading a little on 
15 the top surface according to the preferred embodiment; 

Figure 15 illustrates in a cross -sectional view 
the synchronization of flow in a ribbed plate, with the 
flow splitting and also filling a vertical rib according to 
the preferred embodiment; 
20 Figure 16 illustrates in a cross -sectional view 

the assignment of identical pressures at nodes to 
synchronize flow fronts at a rib according to the preferred 
embodiment; and 

Figure 17 is a view similar to figure 3 generated 
25 according to the preferred eaibodiment of the present 
invention. 

A method for simulating fluid flow within a three 
dimensional object according to the present invention 
30 involves three main steps: 

1) Preparation of geometry, 

2) Selection of injection points, and 

3) Analysis. 

Before these steps are effected, however, the izqput for the 
35 simulation is prepared* Mesh generators for producing 

finite element meshes are generally an expensive addition 
to the base CAD system. Conseguently many CAD 
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installations have no finite element analysis (FEA) mesh 
generation capability. In order to facilitate the 
iniplementation of the invention, the method has been 
designed to avoid the need for a mesh generator. This is 
5 done by using the type of mesh available for 

stereolithography . Stereolithography, a prototyping 
technology that is widely employed in industry, is used to 
produce prototype components for assembly or evaluation. 
It requires a stereolithography apparatus and takes as 

10 input a computer file known as a "SIiA" file. The ability 
to output SLA files is common in CAD systems. However, the 
type of triangulation created is not good enough for FSA. 
COiaracteristic of the SLA file is the inclusion of several 
very long, narrow triangles of high aspect ratio as shown 

15 in figure 4. 

A method in accordance with the present invention includes, 
therefore, a remeshing algorithm (or remesher) that uses as 
input the triangles from a stereolithography sla (or 
20 similar) file and remeshes to produce a mesh that may be 
used in analysis. 

There are several steps involved in the input of a 
stereolithography mesh and its refinement into elements for 
25 analysis. These steps are: 

• Input of mesh, 

• Checking of mesh, 

• Classification of mesh into surfaces, and 

• Remeshing. 

30 Each of these is described below. 

Input of Hesh 

In the preferred embodiment, the input mesh is of the form 
used for stereolithography. Alternatively, however, the 
35 input could be. IGBS surfaces (^^ich would then be meshed 
internally), a surface mesh of quadrilaterals or a mix of 
Quadrilaterals and triangles, tetrahedral or hexahedral 
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meshes • 

Checking of mesh 

The input mesh is checked to ensure that it is closed and 
5 orientable. The mesh is consistently oriented so that a 
unit normal to each element points outward from the 
internal volume defining the solid region that represents 
the object or part being modelled. The volume of this 
region and the surface area of the mesh are calculated ^t**^ 
10 stored. The number of parts defined by the input mesh is 
also determined. 

Classification of mesh into surfaces 

The following properties for each element in the mesh are 
15 noted or calculated, and stored: 

• Area, 

• Normal at the element centroid, 

• Edge lengths of the element, 

• Internal euigles at the vertices of the elements, 
20 • Adjacent elements (the elements attached to each 

side of the element under consideration) , azid 

• Bending angle between adjacent elements (the 
angle between normals of adjacent elements) . 

Nodal properties are then noted, calculated and stored. 
25 These are: ^ 

• Measure of minimum curvature (calculated by 
looking at the elements attached to a node and noting the 
angle between each adjacent pair of elements, whereby the 
measure of minimum curvature is then the smallest such 

30 angle between attached pairs of elements) > 

• Connecting elements (the set of elements 
connected to the node) , and 

• Number of edges connected to the node. 

Surface edges are then calculated using a value of bending 
35 angle to group elements into surfaces. Surfaces so formed 
are then classified into 

• planar surfaces, and 
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• low curvature surfaces. 

Further classification is then done to produce high 
curvature surfaces and edges. The following properties of 
edges are noted or calculated and stored: 
5 • Length, 

• Bending angle, 

• Direction of bending (in or out), and 

• Adjacent Blemants. 

Surface loops are then created. These are defined to be 
10 the oriented edges of surfaces. It should be noted that a 

surface with a hole in it will have two associated loops - 

one for the outer edges and one describing the interior 

hole. Loop properties are then noted or calculated 

stored: 
15 • Length, and 

Bdges connected to loop. 

The following |Burface properties are noted or calculated 

and stored: 

• Perimeter, 
20 • Area, 

• Nodes in surface, 

• Blements in surface, 

• Bdges , and 

• Minimum measure of curvature. 

25 . ' 
Remeshing 

With the surfaces classified, the mesh associated with each 
surface is refined using a bisection algorithm which 
bisects the longest side of a triangle and creates extra 

30 elements. This is illustrated in figures 5a-f . Figure 5a 
shows the initial mesh. The mesh is refined by defining a 
node at the mid point of the longest element side and 
extending lines to vertices to define extra triangles, in 
figures 5b to 5f , the midpoint of the longest side (denoted 

35 by o) and the dotted lines extending from this point define 
the new elements. Remeshing continues until the elements 
satisfy a criterion on size.. 
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The result of the remeshing algorithm when applied to the 
meeh shorn In figure 4 Is showa In figure 6. In this 
regard the ideal triangle is equilateral: it may be seen 
5 that the triangles in figure 6 are far xoore uniform 
closer to this ideal. 

The method will thus accept as Input stereolithography 
files (usually generated with a ".stl" filename extezision) 
10 and other files that are similar such as "render" files 
(produced by "Pro-Bngineer** with a ".sip" extension). 

In practice any file format that describes a mesh (of, for 
example, a lattice or quadrilaterals) covering the outer 
15 surfaces of the solid region may be used as ix^ut to the 

Invention. The remesher is always used to improve the mesh 
quality and produce the triangular mesh for subsequent 
processing . 

20 The first step is then the preparation of the geometry, in 
which the geometry of the three dimensional CAD solid model 
is processed by conventiozial techniques to provide a 
surface mesh covering the outer surfaces of a solid body, 
rather than a midplane mesh as required by the conventional 

25 simulation approach. This distinction Is depicted in 

figures 7a and 7b. Figure 7a shows the midplane mesh for a 
simple plate required by the conventional technique, 
whereas figure 7b shows the mesh raployed according to the 
present invention for the same object. Further, the flow 

30 front must be synchronized on both sides of the surface in 
which the material flows. If material is injected at the 
centre of the plate, then, in the conventional case, the 
niaterial flows to the boundary of the part or object and 
stops without any difficulty. For the surface mesh on the 

35 solid model (figure 7b) the material will flow to the outer 
edge, across the edge and then under the tpp surface. This 
is clearly not Indicative of the physical phenomenon of 
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flow in a plate and so, according to the invention, there 
are imposed some constraints on the flow when using the 
surface mesh. These constraints are that the injection 
point must link to all surfaces from which flow will 
5 emanate and the flow front must be synchronized along the 
top and bottom meshes of the surface mesh. Preferably this 
link is established from the commencement of the flow 
analysis, but a successful (if somewhat less accurate) 
analysis of the flow can still be performed is this linking 

10 commences subsequently. To perform these tasks the surface 
mesh must be categorized into surfaces which are then 
further classified to facilitate injection node selection 
and synchronization. Thus# according to the present 
invention, the surface mesh is then grouped into individual 

15 surfaces and classified. The categories of surfaces are: 

• Matched, 

• Unmatched/ and 

• Edges. 

Hatched surfaces are defined to be those that have a 
20 relationship with another surface such that a notion of 
thickness between them can be sensibly defined. Fisiure 8 
illustrates this idea. Surfaces abed and efgh may clearly 
be associated so that the thickness t between them is 
defined. Surfaces abfe, bcgf , cdhg and daeh are edge 
25 surfaces and unmatched. 

Figure 9 further clarifies these concepts. This figure 
shows a cross-section through a filleted rib. Here ab, cd 
and gh are on edge surfaces. Lines aj and ed are matched 

30 to be, fg and hi . are matched and the curved sections ij and 
ef are unmatched. It should be noted that it is not 
possible to sensibly define a thickness of unmatched 
surfaces in the manner described above. As a fixial 
example, figure 10 illustrates a tapered rib cross-section. 

35 Here ab, cd and gh are edges as in figure 9. Lines aj and 
ed are matched to line be, and the curved lines ef ari/^ ij 
are unmatched as in figure 9. However fg and hi are still 
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considered matched despite the taper. Clearly, if the 
taper is extreme it is possible that the surfaces forming 
the tapered ribs may not be matched. 

5 When surfaces have been classified it is possible to assign 
to elements on matched surfaces a thickness equal to the 
distance between the surfaces. Elements on edge surfaces 
are given the thickness of the matched surfaces to which 
they are attached. Finally elements on unmatched surfaces 
10 are given the average thickness of surrounding elements on 
matched surfaces. 

According to the present invention, the mesh on each side 
of the object is analyzed. A key aspect in this is 

15 obtaining a similar flow front on each side of the matched 
surfaces. As the entry point for flow is the injection 
node and is selected by the user, the method ensures that, 
regardless of where the user chooses to inject plastic, all 
matched surfaces near the point are tied to the injection 

20 node. "Tied" implies that the same pressure and 

temperature of melt at the injection node is given to all 
other tied nodes. 

Figure 11 shows a cross-section of a rectangular plate. 

25 Point A is where the user chooses to inject plastic, point 
AA is the corresponding point on the other side of the 
surface. Point AA is automatically determined in 
accordance with the invention using the information from 
the classification of surfaces described above. Injection 

30 at this point ensures that the flow is matched on top and 
bottom surfaces. Similarly, point B is the point chosen by 
the user and BB is the corresponding point determined by 
the method in accordance with the invention to ensure that 
flow is matched along top and bottom surfaces. 

35 

The selection of injection node(s) can be quite complex. 
This is illustrated in figure 12, which shows an object 
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composed of intersecting plates. The small arrow near the 
point of intersection of the plates indicates the user 
selected injection point. Figure 13 is an enlarged view of 
the injection area and shows the actual injection points 
5 required to allow the flow to emanate correctly from the 
selected point. In this case eight injection points are 
required one of v^ich is supplied by the user and the 
other seven are determined by the invention. 

10 This method identifies which points need to be linked to 
the selected injection point with the concept of matching. 
For a given injection point specified by the user, a set 
ccmsisting of all surfaces which are matched to all 
surfaces coxmected to the injection point is defined. The 

15 following steps are then performed: 

1* For each surface in this set, the point closest 

to the selected injection point is noted. These closest 
points form a set of points that are potentially to be 
linked to the user selected injection point. 

20 2. Each member of the set of potential injection 

points is checked to see if it is already an injection 
point specified by the user, if it has already been so 
determined, it is omitted from the set . 

3. For each potential injection point remaining, a 

25 set consisting of all surfaces which are matched to all 
surfaces connected to the potential injection point is 
defined. 

These steps are repeated until the set of potential points 
30 is unchanged. The set then forms the set of points to be 
linked to the user specified injection point. 

This procedure may be readily understood for the single 
geometry of two crossed plates, for \rtiich the user 
35 specifies injection at the intersection of the plates, in 
the event that the surfaces connected to the specified 
injection point (s) are not matched to any other surfaces. 
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no attenqpt is made to link other points to the injection 
point (s). Instead, flow is allowed to esianate until it 
reaches points attached to surfaces that are matched to 
opposing surfaces. The flow fronts are then synchronised 
5 on the opposing surfaces by assigning links between points 
on the opposing surfaces in a manner similar to that 
described above for linking points to the Injection point. 

Analysis uses a Hele-Shaw type solver. The injection flow 
10 rate is approximately double that for the volume of the 

object as each element in the surface mesh has a thickness 
equal to the gap thickness and so the mesh volume is 
approximately double that of the object. The time step for 
flow front advancement is determined by the geometric 
15 connectivity of the control volwaes attached to each node. 

The solver has the ability to synchronize flow on matched 
surfaces during analysis. Figure 14 shows a cross-section 
of a plate where it is assxmed that the flow front is 

20 advancing from the left and is leading a little on the top 
surface. By definition each node may be empty, partially 
filled or filled. On each time step, Aether all the nodes 
belonging to an element are filled is checked. The instant 
this occurs the element is defined to be the current last 

25 element to fill and a line is drawn from its centrold to 
its matched element (the matched element is known from the 
surface classification described above) defining the check 
point. This "check line** is indicated in figure 14. The 
nearest node to the check point is then assigned the same 

30 pressure as the last node to fill, as iiidicated in figure 
14. Flow syxichronlzatlon is done for all matched surfaces 
and all elements near the flow front (on top and bottom 
surfaces) are checked at each time step. 

35 Synchronization is not performed on unmatched surfaces. 

An example highlighting the role of synchronization is 
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showa in figtire 15, v^ich shows the cross-section of a 
ribbed plate. The flow will split physically and fill the 
vertical rib as well as continuing to the right* Owing to 
the use of the surface aesh it is necessary to impose 
5 synchronization or the flow will continue up through node F 
to the top of the rib and then down to S and C. This is 
performed by using matching. Sizxce the element coxnprising 
nodes D and F and the element comprising nodes C and B will 
be matched (as they are associated with matched surfaces) 

10 it is possible to assign the same pressures at nodes to 

synchronize the fronts. This is illustrated in Figure 16. 
Here a line from the centroid of the element comprising 
nodes D and F is constructed to its matched element 
(comprising nodes C and B) to define the check point. As 

15 the nearest node to the check point is C, it is given the 
same pressure as node D for all subsequent time. Flow now 
emanates from node C as shown and so the flow goes up the 
rib on both sides as required. When establishing links for 
synchronization, links are only made to empty nodes to 

20 which are attached only empty nodes. 

A sample result of the analysis according to the present 
invention is shown in figure 17. This figure corresponds 
to figure 3, which shows the results of the coiqparable, 

25 conventional analysis of the same solid model shown in 
figure 1. As the present invention uses a surface mesh 
rather than a midplane mesh, the results are displayed on 
the solid model (figure 17) rather than on the midplane 
representation as in conventional analysis (figure 3). 

30 Apart from being more natural for the operator, the display 
is visually more meaningful. 

Fundamental quantities calculated by the invention, 
therefore, include: 
35 • the position of the melt front at any time (known 

as "filling contours"), 

• the pressure at the injection node and the 
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pressure distribution throughout the plastic at any time 
during filling or packing phases (known as a ''pressure 
distribution'*}, and 

• the teaic>erature of the plastic melt at any point 
5 within the part at any time during filling or packing 
phases (known as a "temperature distribution"). 
These may be displayed directly or processed to provide 
derived information to the user of the program, such as 
distributions of shear rate, shesu: stress, velocity, 

10 viscosity, direction of flow and orientation of reinforcing 
material. On the basis of this derived Information the 
fundam e ntal calculated quantity data, changes to the 
component geometry, position of injection locations 
(gates), processing conditions used by the injection 

15 molding machine or material for molding may be evaluated 
for their efficacy in improving the quality or 
manufacturability of the part. This is generally an 
Iterative process in \^ich an analysis is performed, 
results are considered and changes made to bring about an 

20 improvement. Another analysis is then performed to ensure 
that the changes have indeed improved the results. 
Frequently changes to the part .geometry are revealed. 
These are made in the CAD system and the revised model is 
subjected to further analysis. It is this aspect which is 

25 assisted by the invention as the need for the designer to 
recreate a new model for analysis every time a change is 
made is removed* 

Thus, an increase in the speed of evaluating of proposed 
30 part shapes is achieved. 

Modifications within the spirit and scope of the invention 
may readily be effected by persons skilled in the art. For 
exan^le, linking may be performed in alternative methods 
35 such as, where multiple Injection points are selected by a 
user, assigning the same flow rate to the linked nodes. It 
is to be xinderstood, therefore, that this invention is not 
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THB CIiAIHS DEFINING TUB INVENTION ARB AS FOLLOWS:. 

1* A method £or alimilating fluid' flow within a three 
dimenBional object having first and second generally 
5 opposed surfaces including: 

matching each element of said first sxirface with 
an element of said second surface between which a 
reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
10 matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 
15 all locations on said first and second surfaces from which 
flow may emanate such that resulting flow fronts along said 
first and second surfaces are synchronized* 

2. A method as claimed in claim 1 wherein said injection 
20 point is first linked to all said locations from 

substantially the commencement of said flow analysis. 

3* A method as claimed in either claim 1 or 2 wherein said 
injection point remains linked to all said locatioxis at 
25 substantially all times in said flow analysis subsequent to 
being first so linked. 

4. A method as claimed in any one of the preceding claims 
wherein said injection point is one of a plurality of 

30 injection points. 

5. A method as claimed in any one of the preceding claims 
wherein said synchronization of said flow fronts is checked 
periodically. 

6. A method as claimed in claim 5 wherein said checking is 
performed at each time step. 
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7 • A method as claimed in any one of the preceding claims 
wherein said £low fronts are synchronized if found not to 
be or no longer to be synchronized. 

5 

8. A method as claimed in any one of the preceding claims 
wherein said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 
surfaces . 

10 

9. A method as claimed in any one of the preceding claims 
wherein any unmatched elements of said first and second 
surfaces, beixig elements that could not be matched, are 
assigned thicknesses being the average of the thicknesses 

15 of adjacent matched elements where such adjacent matched 
el emen ts exist, or of adjacent unmatched elements where 
such adjacent matched elements do not exist ^ry ^ said 
adjacent uxunatc^hed elements have been assigned thicknesses* 

20 10. A method as claimed in claim 9 wherein each element of 
an edge surface, being a surface between said first anH 
second surfaces, and adjacent to either of said first or 
second surface is assigned a thickness proportional to the 
thic3uiess of the element of said first or second surface to 

25 which said element of said edge surface is adjacent. 

11. A method as claimed in claim 10 wherein each said 
element of an edge surface is assigned a thickness between 
0.5 and 1.5 times said thickness of the element of said 

30 first or second surface to which said element of said edge 
surface is adjacent. 

12. A method as claimed in claim 11 wherein each said 
element of an edge surface is assigned a thickness between 

35 0.7 and 0.9 times said thickness of the element of said 

first or second surface to which said element of said edge 
surface is adjacent. 
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13. A metliod as claimed in claim 12 wherein each said 
element of an edge surface is assigned a thickness 
approximately 0.75 times said thiclmess of the element of 

5 said first or second surface to which said element of said 
edge surface is adjacent* 

14 . A method as claimed in claim 10 wherein each element 
of an edge surface not adjacent to said first or second 

10 surface is assigned a thickness being the average of the 
thicknesses of adjacent elements of said edge surface that 
have been assigned thicknesses* 

15. A method as claimed in any one of the preceding claims 
15 wherein flow is simulated at a rate directly proportional 

to a desired flow rate for the object. 

16* A method as claimed in claim 15 wherein said rate is 
proportional to the ratio of computational domain volume of 
20 said object to real volume of said object* 

17* A method as claimed in claim 16 wherein said rate is 
substantially equal to the ratio of said computational 
domain volume to said real volume* 

25 

18. A method as claimed in any one of the preceding claims 
wherein said method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 
30 respectively, wherein said elements are interstices of said 
first and secoxid meshes or lattices. 

19* A method as claimed in any one of the preceding claims 
wherein said elements are triangular* 

35 

20. A method as claimed in any oxie of claims 1 to 18 
wherein said elements are quadrilateral. 
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21. A method as claimed in claim 19 wherein said elements 
are substantially equilateral. 

22. A method as claimed in claim 18 wherein said method 
includes creating said first and second representations. 

23. A method as claimed in either claim 18 or 22 wherein 
said method includes creating improved representations o£ 
said first and second surfaces, whereby said elements are 
elements of said improved representations and said method 
is performed with said improved representations. 



10 



15 



20 



25 



24. A methoid as claimed in claim 18 wherein said first and 
second representations are, or are a part of, a 
representation or representations for stereolithograpl^ of 
said object. 

25. A method as claimed in any one of the preceding claims 
wherein said method is performed by a computer running a 
cwnputer program encoding said method for simulating fluid 
flow. 

26. A method as claimed in any one of the preceding claims 
wherein said method includes corrections for non- isothermal 
temperature fields and/or non-Kewtonian fluids. 



27. A method for simulating fluid flow within a three 
dimensional object having first and second generally 
30 opposed surfaces including: 

providing or creating first and second 
representations of said first and second surfaces 
respectively, 

creating first and second iaiproved 
35 representations from said first and second representations 
respectively, 

matching each element of said first improved 
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representation of said first surface with an element of 
said second improved representation of said second surface 
between which a reasonable thickness may be defined, 
wherein matched elements of said first improved 
5 representation constitute a first set of matched elements 
and matched elements of said second improved representation 
constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using each set of said 
10 matched elements, whereby said injection point is linked to 
all locations on said first and second improved 
representations from which flow may emanate such that 
resulting flow fronts along said first and second improved 
representations are synchronized. 

15 

28. A method as claimed in claim 27 wherein said first and 
second representations are, or are a part of , a 
representation or representations for stereolithography of 
said object. 

20 

29. A method as claimed in either claim 27 or 28 wherein 
said first and second improved representations conprise 
small equilateral triangular elements. 

25 30. A method for simulating fluid flow within a three 
dimensional object having first and second generally 
opposed surfaces including: 

matching each element of said first surface with 
an element of said second surface between which a 

30 reasonable thickness may be defined, wherein matched 

elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements and 
said elements are substantially equilateral triangles, 

35 specifying a fluid injection point, 

performing a flow analysis using each set of said 
matched elements, whereby said injection point is linked to 
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all locations on said first and second surfaces from \^icli 
flow may exaanate such that resulting flow fronts along said 
first and second surfaces are synchronized, wherein said 
first and second representations are, or are a part of, a 
5 representation or representations for stereolithography of 
said object. 

31. A method as claimed in claim 30 wherein said injection 
point is one of a plurality of injection points. 

10 

32. A computing device provided with or running a computer 
program encoding a method for simulating fluid flow as 
claimed in any one of the preceding claims. 

15 33. A computer storage medium provided with a coiqputer 
program embodying a method for simulating fluid flow as 
claimed in any one of claims 1 to 31. 

34. A method for modelling a three dimensional object 
20 including: 

specifying first and secozxd generally opposed 
surfaces of said object, 

forming first and second representations of said 
first axui second surfaces respectively, wherein said first 
25 and second representations each comprise a plurality of 
elements, 

matching pairs of elements of said first and 
second surfaces between which a reasonable thickness may be 
defined. 

30 

35. A method as claimed in claim 34 wherein said first and 
second representations comprise first and second meshes or 
lattices respectively, wherein said elements are 
interstices of said first and second meshes or lattices. 

35 

36. A method as claimed in either claim 34 or 35 wherein 
said elements are triangular. 
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37* A method as claimed in claim 3S wherein said elements 
are substantially equilateral. 

5 38. A method as claimed in claim 34 wherein said elements 
are quadrilateral. 

39* A method as claimed in any one of claims 34 to 38 
wherein each element of each of said matched pairs of 
10 elements is assigned respectively said thickness. 

40. A method as claimed in claim 39 wherein unmatched 
elements of said first and second surfaces are assigned 
thicknesses being the average of the thicScnesses of 

15 surrounding, matched elements of said first and secozxd 
surfaces. 

41. A method as claimed in claim 40 wherein any unmatched 
elements of said first and second surfaces, being elements 

20 that could not be matched, are assigned thicknesses being 
the average of the thicknesses of adjacent matched elements 
where such adjacent matched elements exist, or of adjacent 
unmatched elements where such adjacent matched elements do 
not exist and said adjacent unmatched elements have been 

25 assigned thicknesses. 

42. A method as claimed in claim 41 wherein each element 
of an edge surface, being a surface between said first at ia 
second surfaces, and adjacent to either of said first or 

30 second surface is assigned a thickness proportional to the 
thickness of the element of said first or second surface to 
which said element of said edge surface is adjacent. 

43. A method as claimed in claim 42 wherein each said 

35 element of an edge surface is assigned a thickness between 
0.5 and 1.5 times said thickness of the element of said 
first or second surface to \dxich said element of said edge 



SUBSTITUTE SHEET (RULE 26) 



wo 98/43179 



PCT/AU98/00130 



- 35 - 

surface ia adjacent. 

44. A method as claimed in claim 43 wherein each said 
element of an edge surface is assigned a thickness between 

5 0.7 and 0.9 times said thickness of the element of said 
first or second surface to which said element of said edge 
surface is adjacent. 

45. A method as claimed in claim 44 wherein each said 
10 element of an edge surface is assigned a thickness 0.75 

times said thickness of the element of said first or second 
surface to which said element of said edge surface is 
adjacent. 

15 46. A method as claimed in claim 42 wherein each element 
of an edge surface not adjacent to said first or second 
surface is assigned a thickness being the average of the 
thicknesses of adjacent elements of said edge surface that 
have been assigned thicknesses. 

20 

47 . A method for simulating fluid flow within a three 
dimensional object having first and second generally 
opposed surfaces including: 

matching each element of said first surface with 
25 an element of said second surface between which a 

reasonable thickness may be defined, wherein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements , 
30 specifying a fluid injection ipoint, 

performing a flow analysis using each set of said 
matched elements, and 

synchronizing flow fronts resulting from said 
flow analysis along said first and second surfaces. 
35 * 

48. A method as claimed in claim 47 wherein said flow 
fronts are synchronized from substantially the commencement 
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o£ said flow analysis. 

49. A method as claimed in claim 47 wherein said flow 
fronts are first synchronized after the commencement of 

5 said flow analysis. 

50. A method as claimed in any one of claims 47 to 49 
wherein said in:j action point remains linked to all said 
locations at substantially all times in said flow analysis 

10 subseguent. to being first so linked.. 

51. A method as claimed in any one of claims 47 to 50 
wherein said injection point is one of a plurality of 
injection points. 

15 

52. A method as claimed in any one of claims 47 to 51 
wherein said synchronization of said flow fronts is checked 
periodically. 

20 53. A method as claimed in claim 52 wherein said checking 
is performed at each time step. 

54. A method as claimed in any one of claims 47 to 53 
wherein said flow fronts are syachronized if found not to 

25 be or no longer to be synchronized. 

55. A method as claimed in any one of claims 47 to 54 
wherein said first and second generally opposed surfaces 
are one of a plurality of pairs of generally opposed 

30 surfaces. 

56. A method as claimed in any one of claims 47 to 55 
wherein any unmatched elements of said first and second 
surfaces, being elements that could not be matched, are 

35 assigned thicknesses being the average of the thicknesses 
of adjacent matched elements where such adjacent matched 
elements exist, or of adjacent unmatched elements where 
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Buch adjacent matched elements do not exist said 
adjacent unmatdxed elements have been assigned thicknesses. 

57. A method as claimed in claim 56 wherein each element 

5 of an edge surface, being a surface between said first and 
second surfaces, and adjacent to either of said first or 
second surface is assigned a thickness proportioxial to the 
thickness of the element of said first or second surface to 
which said element of said edge surface is adjacent. 

10 

58. A method as claimed in claim 57 wherein each said 
element of an edge surface is assigned a thickness between 
0.5 and 1.5 times said thickness of the element of said 
first or second surface to which said element of said edge 

15 surface is adjacent. 

59. A xoethod as claimed in claim 5B wherein each said 
element of an edge surface is assigned a thickness between 
0.7 and 0.9 times said thickness of the element of said 

20 first or second surface to which said element of said edge 
surface is adjacent. 

60. A method as claimed in claim 59 wherein preferably 
each said element of an edge surface is assigned a 

25 thickness approximately 0.75 times said thickness of the 
element of said first or second surface to which said 
element of said edge surface is adjacent. 

61. A method as claimed in claim 60 wherein each element 
30 of an edge surface not adjacent to said first or second 

surface is assigned a thickness being the average of the 
thicknesses of adjacent elements of said edge surface that 
have been assigned thicknesses. . 

35 62. A method as claimed in any one of claims 47 to 61 

wherein flow is simulated at a rate directly proportional, 
to a desired flow rate for the object. 
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63. A metlxod as claimed in claim 62 \^erein said rate is 
proportional to the ratio of conqputational domain volume of 
said object to real volume of said object. 

.5 

64. A method as claimed in claim 63 therein said rate is 
substantially equal to the ratio of said computational 
domain volume to said real volume. 

10 65. A method as claimed in any one of claims 47 to 64 
wherein said method is performed with first and second 
representations of said first and second surfaces 
respectively comprising first and second meshes or lattices 
respectively, herein said elements are interstices of said 

15 first and second meshes or lattices. 

66. A method as claimed in any one of claims 47 to 65 
wherein said elements are triangular. 

20 67. A method as claimed in any one of claims 47 to 65 
wherein said elements are quadrilateral. 

68. A method as claimed in claim 66 wherein said elements 
are substantially equilateral. 

25 

69. A method as claimed in claim 65 wherein said method 
includes creating said first and second representations. 

70. A method as claimed in claim 65 wherein said method 
30 includes creating improved representations of said first 

and second surfaces, whereby said elements are elements of 
said improved representations and said method is performed 
with said improved representations. 

35 71. A method as claimed in claim 65 idierein said first and 
second representations are, or are a peirt of, a 
representation or representations for stereolithography of 
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said Object. 

72. A aethod aa daisied in any one of claims 47 to 71 
<«ierein said atethod includes corrections for non- isothermal 

5 teaperature fields and/or non-Newtonian fluids. 

73. A method as claimed in any one of claims 47 to 72 
wherein said method is performed by a coavuter running a 
conqputer program encoding said method for simulating fluid 

10 flow. 

74. A method for modelling three dimensional objects by 
finite element analysis including: 

modifying a representation for stereolitixography 
15 of said object to a representation suitable for finite 
element analysis. 

75. A method as claimed in claim 74 wherein said method 
includes creating improved representations of said first 
and second surfaces, whereby said elements are elements of 
said improved representations and said method is performed 
with said improved representations. 

76. A method for simulating fluid flow within a three 
25 dimensional object having first and second generally 

opposed surfaces including: 

matching each element of said first surface with 
an element of said second surf ace between ^^ich a 
reasonable thickness may be defined, idierein matched 
elements of said first surface constitute a first set of 
matched elements and matched elements of said second 
surface constitute a second set of matched elements, 
specifying a fluid injection point, 
performing a flow analysis using said first set 
35 of matched elements, 

adapting said flow analysis to said second set of 
matched elements, and 



20 



30 
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synchronizing flow fronts resulting from said 
flow analysis and said adaptation of said flow analysis 
along said first and second surfaces. 

5 77. A method as claimed in claim 76 wherein said "method is 
performed with first and second representations o£ said 
first and second surfaces respectively conprising first n n d 
second xoeshes or lattices respectively, wherein said 
elements are interstices of said first and second meshes or 
10 lattices . 

78. A method as claimed in either claim 76 or 77 wherein 
said elements are triangular. 

15 79. A method as claimed in either claim 76 or 77 wherein 
said elements axe quadrilateral. 

80. A method as claimed in claim 78 wherein said elements 
are substantially equilateral. 

20 

81. A method as claimed in any one of claims 76 to 80 
wherein said method includes creating said first and second 
representations . 

25 82. A method as claimed in any one of claims 76 to 81 
wherein said method includes creating improved 
representations of said first and second surfaces, whereby 
said elements are elements of said improved representations 
and said method is performed with said iinproved 

30 representations. 

83. A method as claimed in any one of claims 1, 27 or 47 
wherein said synchronization comprises matching pressure 
axLd temperature. 
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Figure 6 
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Figure 12 
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Figure 13 
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